Design of Metamaterial Surfaces with Broad-band Absorbance 
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A simple design paradigm for making broad-band ultra-thin plasmonic absorbers is introduced. 
The absorber's unit cell is composed of sub-units of various sizes, resulting in nearly 100% absorbance 
at multiple adjacent frequencies and high absorbance over a broad frequency range. A simple 
theoretical model for designing broad-band absorbers is presented. It uses a single-resonance model 
to describe the optical response of each sub-unit and employs the series circuit model to predict 
the overall response. Validity of the circuit model relies on short propagation lengths of the surface 
plasmons. 
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Plasmonic metamaterials (MMs) are subwavelength 
metallic nanostructures designed to provide extraordi- 
nary electromagnetic properties in specified frequency 
ranges [1, 2]. The combination of the plasmonic response 
of conduction electrons and appropriate (but not neces- 
sarily geometrically complicated) design of the structures 
results in subwavelength feature sizes and strong near- 
field enhancement. These properties pave the way to var- 
ious applications, such as surface enhanced Raman scat- 
tering [3, 4], light trapping [5], bio- molecule sensing [6], 
and plasmonic absorbers [7-12]. For the reasons going 
beyond the simplicity of fabrication, many applications 
utilize MM surfaces instead of MMs. For example, MM 
surfaces provide solid/liquid interfaces for bio-sensing, or 
can be used as ultra-thin coatings to provide spectrally- 
selective optical properties (e.g, absorption). 

Simple designs of narrow-band MM absorbers 
(NBMA) comprised of an array of plasmonic 
strips/patches separated by a thin dielectric spacer 
from a ground plate, as illustrated in Fig. 1(a), have 
been recently proposed and experimentally imple- 
mented [8, 10-12]. "Perfect" (100% at the peak) light 
absorption by NBMAs is enabled by the resonant 
excitation of a surface plasmon polariton (SPP) whose 
frequency is controlled by the structural dimensions. 
Using such NBMAs as building blocks for a broad-band 
metamaterial absorber (BBMA) appears straightfor- 
ward: combining several such blocks (i — 1, ...N) with 
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In reality, designing BBMAs using this approach [13] 
has been challenging. For example, the individual 
absorption peaks of the constitutive NBMAs can each 
yield "perfect" absorption, but at the expense of large 
spectral separation. Alternatively, constitutive NBMA's 
with spectrally-close absorption peaks can fail to provide 
"perfect" absorption at every peak. Systematic approach 
to designing NBMAs and combining them into a unit 
cell of a BBMA is needed. In this Letter, we provide 
such an approach based on a simple impedance model. 
We show that the surface impedance z of a BBMA can 
be obtained from the impedances Zi of the constitutive 
NBMAs using a simple additive formula (see Eq. (3)), 
which is valid as long as the SPPs' propagation lengths 



are much smaller than the separation between NMBAs. 
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FIG. 1: (a) Schematic of the Narrow Band Metamaterial 
Absorber (NBMA). (b) Dependence of radiative (wj e ) and 
resistive (wi a ) decay rates on G. Remaining dimensions: 
L = 350 nm, W = 250 nm, and D = 20 nm. Inset: resonant 
wavelength vs G. (c) Impedances of NBMAs for 3 values of 
G. Solid lines: from driven simulations, dashed lines: from 
the single oscillator model, (d) Absorbance of the NBMAs 
with same values of G as in (c). 

We start by introducing a simple single-resonance 
model of an NMBA based on the periodic structure 
shown in Fig. 1(a). The structure coupled to the inci- 
dent radiation is described within the framework of a 
single resonator model [14]. An open resonator is de- 
fined by its natural frequency, uq, and finite lifetimes, 
t = 1/uJio and -r e = l/wie, respectively determined by 
the Ohmic and radiative losses. The interaction between 
the resonator and the incident field of amplitude Ej is 
described by the following equations: 



—a = -iuj a - (wj + CJ ie ) a + \J2 uj ie Ej, 
at 



E R = -Ei + V2 Lo ie a, 



(1) 



where a and En are the amplitudes of the resonator and 
the reflected field, respectively. The specific geometry of 
the resonator that includes the ground plate is manifested 
in the expression for Er, which is a superposition of the 
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resonator and ground plate reflection. Within this simple 
single-resonance model, the surface impedance defined as 
z = (1 + r)/(l — r), where r = En/Ej, is given by: 

(2) 



—i(u - cj ) + 

and all optical properties of the NBMA are determined 
by three parameters: uq, uji e , and Wj D . These parame- 
ters depend on the structural dimensions and can be ob- 
tained by calculating the "leaky" SPP eigenmodes [12] 
using eigenvalue simulations (i.e., no incident field). By 
tuning the geometric parameters of the NBMA, one can 
engineer u>q, Ui e , and L0i O and achieve z = 1 (and, there- 
fore, perfect absorption) whenever the critical coupling 
condition LOi e = uji is satisfied. 

The validity of the single-resonance model was ver- 
ified for several NBMAs as shown in Fig. 1. Plas- 
monic strips and ground plate are assumed to be made 
of gold (described by the Drude permittivity e = 1 — 
u)p/ [u>(u> + icj c )], where u> p — 1.32 x 10 w rad/s and u> c — 
1.2 x 10 u rad/s) and the spacer is assumed to have a 
non-dispersive dielectric permittivity = 2.25. The ac- 
curacy of the z(bS) obtained from Eq. (2) was verified 
by comparing it to the impedance obtained from the re- 
flection coefficients calculated at various frequencies w 
of the incident radiation. Excellent agreement shown 
in Fig. 1(c) validates the single-frequency model which 
is highly computationally efficient because only a single 
eigenvalue/ eigenmode simulation is required [12] to de- 
termine z(lu) for all frequencies. 

The tunability of the NBMA's frequency characteris- 
tics (ujq, u ie , uj io ) with respect to the spacer thickness G 
is illustrated in Fig. 1(b). Note from Fig. 1(d) that con- 
siderable peak absorbances (in excess of 90%) is achieved 
for a wide range of spacer thicknesses lOnm < G < 20nm. 
The smooth dependence of the frequency characteristics 
on G shown in Fig. 1(b) (as well as on the other struc- 
ture parameters) enables rapid design of the appropriate 
parameter set (W, L, G, D) of the NBMA that provides 
the targeted (w , w ie , uj io ). 

Two important features of the NBMAs described in 
Fig. 1 make them promising building blocks for de- 
signing a BBMA. First, their period is strongly sub- 
wavelength: Li <C Aq 1 ' = 2ttc/luq \ That implies that 
a sub-wavelength (L — £\ Z, < A) unit cell of a BBMA 
can accommodate a large number of NBMA-based sub- 
units satisfying Li <C L < A. Second, the propagation 
length igpp of the SPP responsible for absorption satisfies 

4pp < This property has been shown [12] to be re- 
lated to the wide-angle absorbance of the NBMA. While 
these two conditions are related, they are not equivalent 
as shown below. Nevertheless, when both conditions are 
satisfied, a simple expression for a the surface impedance 
of a BBMA with a macro-cell composed of N NBMA's 
can be obtained: 



N 



where (Li) is the impedance (period) of each periodic 
constitutive NBMA. Derivation of Eq. (3), based on the 
series combination model [15], is illustrated in Fig. 2(a), 
where the external field represents a voltage drop and 
each NBMA resonator is treated as a circuit element. 
The key assumption is that each component of the cir- 
cuit functions independently from each other, i.e., there 
is little mutual capacitive or inductive coupling between 
adjacent sub-units. This condition is satisfied owing to 
the second feature of NBMAs. Below we verify Eq. (3) 
for two cases: (i) strongly (Fig. 2(b,c)) and (ii) weakly 
(Fig. 2(d,e)) confined SPPs of the individual sub-units. 
In each case the unit cell of a candidate BBMA consists 
of two sub-units with dimensions stated in Fig. 2. 
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FIG. 2: (a) Illustration of the circuit model used to derive 
Eq. (3). (b,d) Impedance z of the resonant surface with a 
macro-cell comprised of two sub-units. Sub-unit dimensions: 
(b) [Li, L a , Wi, Wa] = [407nm, 383nm, 268nm, 227nm] and 
[D,G] = [17nm, 28nm]; (d) [L 1; L 2 , W u W 2 ] = [432nm, 371nm, 
409nm, 356nm] and [D,G] = [13nm, 15nm]. Solid lines in (b,d): 
z from direct simulations; dashed lines: z from Eq. (3), with 
individual impedances z\^, plotted in (c,e), respectively. 

In case (i), excellent agreement between BBMA's 
impedance calculated from Eq. (3) and from direct nu- 
merical simulations is shown in Fig. 2(b). Each sub-unit, 
representing a resonant high- impedance (zi,2 ~ 2) sur- 
face, is independent of the other inside the macro-cell 
because ?spp ~ 360nm < L\^.- The impedance of the 
combined candidate BBMA has two spectral peaks with 
values z w 1 corresponding to 100% absorption. On the 
other hand, Eq. (3) clearly fails in case (ii) as shown in 
Fig. 2(d). That is because ^spp ~ 705nm > L,-, and 
the two sub-units strongly interact inside the macro-cell, 
thereby invalidating the series impedance model. Note 
that to achieve such long propagation length of the SPP, 
the electron collision frequency, u) c , in Au has been arti- 
ficially reduced by a factor of four. 

We now present two examples of BBMAs constructed 
by engineering their unit cells from several (two in 
Fig. 3(a,c) and three in Fig. 3(b,d)) NBMA sub-units. 
We use the example of two sub-units inside the unit 
cell of a BBMA to describe the computational proce- 
dure for obtaining an absorber with two 100% absorp- 
tion peaks at the prescribed frequencies uj\ and u>2- That 
corresponds to z(u>{) — ziyi?) — 1, where z(u>) is calcu- 
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FIG. 3: Absorption spectrum (top row) and surface 
impedance (bottom row) of the designed Broad Band Meta- 
material Absorbers (BBMA) with two (left column) and three 
(right column) sub-units per unit cell. Sub-units dimensions: 
[Li, L 2 ,W ly W 2 ] = [448nm, 348nm, 249nm, 226nm] in (a,c); 
[Li, L 2 , L 3 ] = [327nm, 280nm, 271nm] and [Wi, W 2 , W 3 ] = 
[224nm, 200nm, 180nm] in (b,d). Metal and spacer thickness 
[D, G] = [17nm, 28nm] is the same for all sub-units. 

latcd from Eq. (3). The resulting four equations (note 
that z is a complex number) must be satisfied for four 
unknown. For simplicity, we assume D\ = D2 and 
Gi = G2 to be fixed. This leaves the complex impedance 
z = z{L\, L21 Wi, W2) as a function of four unknowns 
that can be numerically solved for to satisfy the above 
four equations. The numerical implementation based on 
the multi-dimensional Newton's method is greatly sim- 
plified by the variables' separation in Eq. (3): the num- 
ber of partial derivatives dz/dLi, dz/dWi is reduced be- 
cause dzi/d(L,W)j = 5ijdzo/d(L,W), where zq is the 
impedance of one single sub-unit. 

The results of such numerical calculations for a BBMA 



with two sub-unit per unit cell are shown in Fig. 3 (left 
column) , where the resulting double-peaked impedance z 
is plotted in Fig. 3(c) and the corresponding frequency- 
broadened absorptivity is plotted in Fig. 3(a). A simi- 
lar numerical procedure was employed to design a triple- 
peaked BBMA illustrated in Fig. 3 (right column; all 
sub-unit parameters in the caption). Note that over 
90% absorptivity is achieved over the spectral range of 
1.2/xm < A < 1.5fim. Note that in both examples pre- 
sented in Fig. 3 the periodicity Lbbma = J2i Li of the 
BBMA is sub-wavelength: Lbbma < A. It is essential to 
satisfy this condition to avoid the emergence of the ad- 
ditional diffraction orders which create additional radia- 
tive loss channels and further complicate the impedance 
model. In practice, this sets a limit on the attainable 
bandwidth. For example, for the strip-based designs pre- 
sented in this Letter, the number of sub-units is limited 
to 3 (or 9 for a patche-based design). A possible solu- 
tion could involve an aperiodic distribution of sub-units. 
Based on negligible interaction between the sub-units, it 
is expected that Eq. (3) should hold. 

In conclusion, we have demonstrated that a broad- 
band metamaterial absorber (BBMA) can be designed 
by constructing a super-lattice structure with each sub- 
lattice absorbing at different frequencies. Based on a 
single oscillator model and a series circuit model, broad- 
band absorbers with multiple perfectly absorbing peaks 
are shown to be attainable, and a procedure for design- 
ing such absorbers is provided. Examples of BBMAs 
comprised of two and three sub-units operating in the 
telecommunications frequency range are presented, and 
an average 90% absorbance in the 1.2 fim < A < 1.5/im 
range is predicted. Such ultra-thin (< lOOnm) BBMAs 
can be used for infrared detection and other applications. 
This work is supported by the Office of Naval Research 
(ONR) grant N00014-10-1-0929. 
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